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Abstract 

Active walker models have recently proved their great value for describing the 
formation of clusters, periodic patterns, and spiral waves as well as the devel- 
opment of rivers, dielectric breakdown patterns, and many other structures. 
It is shown that they also allow to simulate the formation of trail systems 
by pedestrians and ants, yielding a better understanding of human and ani- 
mal behavior. A comparison with empirical material shows a good agreement 
between model and reality. 

Our trail formation model includes an equation of motion, an equation for 
environmental changes, and an orientation relation. It contains some model 
functions, which are specified according to the characteristics of the considered 
animals or pedestrians. Not only the kind of environmental changes differs: 
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Whereas pedestrians leave footprints on the ground, ants produce chemical 
markings for their orientation. Nevertheless, it is more important that pedes- 
trians steer towards a certain destination, while ants usually find their food 
sources by chance, i.e. they reach their destination in a stochastic way. As a 
consequence, the typical structure of the evolving trail systems depends on the 
respective species. Some ant species produce a dendritic trail system, whereas 
pedestrians generate a minimal detour system. 

The trail formation model can be used as a tool for the optimization 
of pedestrian facilities: It allows urban planners to design convenient way 
systems which actually meet the route choice habits of pedestrians. 



Typeset using REVT^ 
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I. INTRODUCTION 
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The emergence of complex behavior in a system consisting of simple, interacting elements 
is among the most fascinating phenomena of our world. Examples can be found in 
almost every field of today's scientific interest, ranging from coherent pattern formation in 
physical and chemical systems to the motion of animal swarms in biology and 

the behavior of social groups [§-[n]|. 

In the life and social sciences, one is usually convinced that the evolution of social sys- 
tems is determined by numerous factors, such as cultural, sociological, economic, political, 
ecological etc. However, in recent years, the development of the interdisciplinary field "sci- 
ence of complexity" has lead to the insight that complex dynamic processes may also result 
from simple interactions, and even social structure formation could be well described within 
a mathematical approach |jlO|-0]. Moreover, at a certain level of abstraction, one can find 
many common features between complex structures in very different fields. 

A recent field of particular interest is the microsimulation of self-organization phenomena 
occuring in traffic systems. This includes the formation of jammed states in freeway or city 



traffic |15-28| as well as the various collective patterns of motion developing in pedestrian 



crowds p8|-|32| like oscillatory changes of the walking direction at narrow passages or round- 
about traffic at crossings. 

In this paper, we draw the attention to the specific collective phenomenon of trail for- 



mation |l33|j3^ , which is widely spread in the world of animals and humans. Regarding their 
shape, duration and extension, trail systems of different animal species and humans differ, 
of course. However, more striking is the question, whether there is a common underlying 
dynamics which allows for a generalized description of the formation and evolution of trail 
systems. 

As our experience tells us, trails are adapted to the requirements of their users. In the 
course of time, frequently used trails become more developed, making them more cittrcLC- 
tive, whereas rarely used trails vanish again. Trails with large detours become optimized 
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by creating shortcuts. New destinations or entry points are connected to an existing trail 
system. These dynamical processes occur basically without any common planning or direct 
communication among the users. Instead, the adaptation process can be understood as a 
self-organization phenomenon, resulting from the non-linear feedback between the users and 



the trails BH]. 



In order to simulate this process, we propose here a particle-based, multi-agent approach 
to structure formation, which belongs to the class of active walker models. Like random walk- 
ers, active walkers are subject to fluctuations and influences of their environment. However, 
they are additionally able to locally change their environment, e.g. by altering an environ- 
mental potential, which in turn influences their further movement and their behavior. In 
particular, changes produced by some walkers can influence other walkers. Hence, the non- 
linear feedback can be interpreted as an indirect interaction between the active walkers via 
environmental changes, which may lead to the self-organization of spatial structures. 

Active walker models have proved their versatility in a variety of applications, such 
as formation of complex structures p6|-|4^, pattern formation in physico-chemical systems 
^3|^6[, aggregation in biological P7| , ^ or urban systems, and generation of directed 



motion [BD|,ET|. The approach provides a quite stable and fast numerical algorithm for sim- 



ulating processes involving large density gradients, and it is applicable also in cases where 
only small particle numbers govern the structure formation. In particular, the active walker 
model is applicable to processes of pattern formation which are intrinsically determined by 
the history of their creation, such as the formation of trail systems, discussed in this paper. 

In Section II, the active walker model for trail formation is formulated in terms of a 
Langevin equation for the movement of the walkers, an equation for environmental changes, 
and a relation describing the orientation of the walkers with respect to existing trails. As 
one application of the model. Section HI describes the formation of trunk trails in ant 
colonies, which are commonly used to exploit food sources. As a second application, in 
Section IV the evolution of pedestrian trail systems is modelled. Both Sections HI and IV 
present a comparison of computational results with real trail systems, indicating a good 
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aggrement between model and empirical facts. In Section IV.A, the equations for pedestrian 
trail systems are scaled to dimensionless equations, in order to demonstrate that the evolving 
trail systems are (apart from the boundary conditions) only determined by two parameters. 
In Section IV. B, a macroscopic formulation of human trail formation is derived from the 
microscopic equations, allowing analytical investigations and an efficient calculation of the 
stationary solution by a self-consistent field method. Our conclusions and an outlook, which 
suggests an application of the model to the optimization of trail systems, are presented in 
Section V. 

II. ACTIVE WALKER MODEL OF TRAIL FORMATION 

In order to introduce our model, we first describe the process of trail formation within 
a general stochastic framework. Hence, in this section the active walkers are not specified 
as pedestrians or animals. Rather, they are considered as arbitrary moving agents, who 
continuously change their environment by leaving markings while moving. These markings 
can, for example, be imagined as damaged vegetation on the ground (as in the case of hoofed 
animals or pedestrians) or as chemical markings (as in the case of ants). 

The spatio-temporal distribution of the existing markings will be described by a ground 
potential Gk{r,t). Trails are characterized by particularly large values of Gk{r,t). The sub- 
script k allows to distinguish different kinds of markings. Due to weathering or chemical 
decay, the markings have a certain life time Tk{r) which characterizes their local durability. 
Therefore, existing trails tend to fade, and the ground potential would exponentially adapt to 
the natural ground conditions G^(r), if the production of markings would be stopped. How- 
ever, the creation of new markings by agent a is described by the term Qai^a, t)S{r — r^), 
where Dirac's delta function S{r — Va) gives only a contribution at the actual position ra{t) 
of the walker. The quantity QaiTa-it) represents the strength of new markings and will 
be specified later. In summary, we obtain the following equation for the spatio-temporal 
evolution of the ground potential: 
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■[Gl{r)-G,ir,t)] 



dGk{r,t) _ 1 



dt Tkir) 

+ 5^Qa(r-a,t)5(r-r„(t)). (1) 

a 

The motion of the active walker a on a two-dimensional surface will be described by the 
following Langevin equation: 

dVait) 



dt 

dVa(t) 



Va{t) (2a) 



-laVait) + /Jt) + v/2£„7aCW- (2b) 



dt 

Eq. (0) considers both deterministic and stochastic influences on the motion of the active 
walker. Va denotes the actual velocity of walker a. 7q, represents some kind of friction coeffi- 
cient. It is given by the relaxation time of velocity adaptation, specified later: 7^ = l/r^. 
The last term describes random variations of the motion in accordance with the fluctuation- 
dissipation theorem. Ea is the intensity of the stochastic force ^^(t), which was assumed to 
be Gaussian white noise: 

iUt)) = , iUmp^it')) = 5.^5,,5{t - t'). (3) 

The a-dependence of Sa takes into account that different walkers could behave more or less 
erratic, dependent on their current situation. 

Finally, the term represents deterministic influences on the motion, such as intentions 
to move into a certain direction with a certain desired velocity, or to keep distance from 
neighboring walkers. According to the social force concept |2^j3^ , /„ is specified as follows: 

^0 

fait) = —ea{ra,Vo„t) 

+ Yl fa(5{ra,v^;rp,vp) . (4) 

Here, f ^ describes the desired velocity and the desired direction of the walker. The term 
f delineates the effect of pair interactions between walkers a and (3 on the motion of 
walker a P5|J5I1|J53[ ] . Since we will focus on cases of rare direct interactions, f^^ can be 
approximately neglected here. Thus, eq. (||b) becomes 
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dVa(t) 

dt 




(5) 



TV 



a 



where the first term reflects an adaptation of the actual walking direction Va/Hi^all to the 
desired walking direction and an acceleration towards the desired velocity with a 
certain relaxation time r^. Assuming that the time Tq is rather short compared to the time 
scale of trail formation (which is characterized by the durability T^), equation can be 
adiabatically eliminated. This leads to the following equation of motion: 



To complete our trail formation model, we must finally specify the orientation relation 



which determines the desired walking direction in dependence of the ground potentials 
Gk{r,t). Since the concrete orientation relation for pedestrians differs from that for ants, it 
will be introduced later on, in the respective sections. However, it is clear that the presence 
of a trail will have an attractive effect, i.e. it will induce an orientation towards it. According 



Therefore, the mechanism of trail formation is based on some kind of agglomeration pro- 
cess, which is delocalized due to the directedness of the walkers' motion. Starting with a 
plain, spatially homogeneous ground, the walkers will move arbitrarily. However, by con- 
tinuously leaving markings, they produce trails which have an attractive effect on nearby 
walkers. Thus, the agents begin to use already existing trails after some time. By this, a 
kind of selection process between trails occurs (cf. ^^): Frequently used trails are reinforced, 
which makes them even more attractive, whereas rarely used trails may vanish again. The 
trails begin to bundle, especially where different trails meet or intersect. Therefore, even 
walkers with different entry points and destinations use and maintain common parts of the 
trail system. 






(7) 



to Eq. (I), this will cause a tendency to approach and to use the trail. 
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As a first example, we want to model the formation of trunk trails, which is a widely ob- 
served phenomenon in ant colonies, such as in the Myrmicinae, Dolichoderinae and Formic- 
inae species, commonly foraging for food from a central nest [^-53|. The trails are used 
to connect the food sources with the nest to allow for a collective exploitation of the food. 
In the case of ants, the markings are chemical signposts, so-called pheromones, which also 
provide the basic orientation for foraging and homing of the animals. However, note that 
not all ants species form trails. There is a variety of very complex foraging patterns in ants, 
such as swarm riding of army ants (e.g. in the species of Eciton and Dorylus) [Q. Therefore, 
we restrict here to cases, in which trunk trail formation of group riding ants is reported. 

Before we present our model, we would like to mention some differences between active 
walkers and ants. The latter are rather complex biological creatures which are capable of 



using additional information (e.g. landmark use) or egocentric navigation for their food 
searching and homing. Moreover, they can store information in an individual memory and 
communicate with nest mates in a very complex manner |^ . 



We will neglect these abilities, in order to show that they are not necessary for trail 
formation. The active walkers in our model merely count on the local information provided 
by the chemical trail, in order to guide themselves. They do not have additional navigation 
or information processing capabilities, and are not subject to long-range attracting forces to 
the food sources or to the nest. Hence, the formation of trunk trails in the following model 



is clearly a self-organizing process, based on the local interactions of the walkers fST 



Trunk trails used for foraging are typically dendritic in form. Each one starts from the 
nest vicinity as a single thick pathway that splits first into branches and then into twigs to 
convey large numbers of ants rapidly into the foraging areas (see Fig. 1). 

In order to distinguish those trails which lead to a food source, the ants, after discovering 
a food source, use another pheromone to mark their trails, which stimulates the recruitment 
of additional ants to follow that trail. In our active walker model, we count on that fact 
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by using two different chemical marlcings: Chemical is used by the active walkers as long 
as they have not reached a food source, i.e. on their way from the nest to the food or 
during search periods. Chemical 1 is only used by active walkers after they have reached 
a food source, i.e. on their way back from food sources to the nest. An internal parameter 
ka = {0, 1} indicates which of these markings is produced by the active walker a. Hence, 
the production term for the ground potential in Eq. ([l|) is defined as follows: 

Qaira, t) = (1 - ka)qo exp[-/3o (t - to)] 

+fcagiexp[-/5i(t-t?)]. (8) 

The first term is relevant for = 0, i.e. when searching a food source, whereas the second 
term contributes only for ka = 1, i.e. after having found some food. Since the capacity 
of producing chemical markings is limited, we have assumed that the quantity of chemical 
produced by a walker after leaving the nest or the food source decreases exponentially in 
time, where Po and Pi are the respective decay parameters, go and qi denote the initial 
production, and tg, are the times, when the walker a has started from the nest or the 
food source, respectively. 

Due to the two chemical markings, we have two different ground potentials Go{r,t) and 
Gi{r,t) here, which provide orientation for the walkers. In the following, we need to specify 
how they influence the motion of the agents, especially their desired directions ea{ra, Va, t). 
At this point, we take into account that the walkers a will not directly be affected in their 
behavior by the ground potentials Gk{r,t) themselves, which reflect the pure existence of 
markings of type k at place r. They will rather be influenced by the perception of their 
environment from their actual positions r^, which will be described by the trail potentials 

V,t{ra, t) = V,ti{G,{r, t)}, r„, v^) . (9) 

For the detection of chemical markings, insects like ants use specific receptors which are 
located at their so-called antennae. Their perception is mainly determined by the angle 2ip 
of perception, which is given by the angle between the antennae (cf. Fig. 2). Therefore, we 
make the assumption 
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dr' / dip' r' 
J 

X Gk [va + r'{ cos{uJa + ^'), sin(u;„ + V^'))) ^) > (10) 

where the angle Ua is given by the current walking direction 

= ^77M = (cos.;.(t),sin.;4t)). (11) 

According to (0), our active walkers integrate over the ground potential between the an- 
tennae of length Ar. Note, however, that the restriction to the angle of perception is not an 



indispensible assumption for the generation of trails fS^]- Thus, it could be neglected in a 
minimal model. Nevertheless, it has been introduced to mimic the biological constitution of 
the ants and to keep close to biology. 

The perception of already existing trails will have an attractive effect ftri'^^ai i^a-, t) to the 
active walkers. This has been defined by the gradients of the trail potentials, 

fUr,v,t) = {1 - k^)VV,l{r,v,t) 

+KVV;^,{r,v,t). (12) 

The above formula takes into account that walkers which move out from the nest to reach 
a food source {ka = 0) orientate by chemical 1, whereas walkers which move back from the 
food {ka = 1) orientate by chemical 0. This implies that initially, in the absense of chemical 
1, the walkers move as random walkers which discover a food source only by chance. 

We complete our model of trunk trail formation by specifying the orientation relation of 
the walkers. Assuming ea{r,v,t) = /°^(r, v, t)/||/^(r, v, the desired walking direction 
ea{r,v,t) points into the direction of the steepest increase of the relevant trail potential 
V^^(r,f,t). However, this formula does not take into account the ants' persistence to keep 



the previous direction of motion ||58|. The latter reduces the probability of changing to the 
opposite walking direction by fluctuations, which would cause the ants to move backwards 
before reaching their goal. Therefore, we modify the above formula to 
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ea{r,v,t) 



(13) 



where Naif.v.t) = \\f%.{r,vj) + e*(t — At)|| is a normalization factor. That means, on a 
ground without markings, the walking direction tends to agree with the one at the previous 
time f — At, but it can change by fluctuations. 

Finally, it is known from ant species that they are able to leave a place where they do 
not find food and increase their mobility to reach out for other areas. Since active walkers 
do not reflect their situation, they stick on their local markings even if they did not find any 
food source. In order to increase the mobility of the active walkers in those cases, we assume 
that every walker has an individual noise intensity £a(t), which is related to the walker's 
spatial diffusion coefficient and should increase continuously, as long as the walker does not 
find a food source: 



is again the starting time of walker a from the nest, Sq is the initial noise level and its 
growth rate. If the noise intensity Sa has reached a critical upper value, the walker a behaves 
more or less as a random walker which does not pay attention to the trail potential. But if 
the walker found some food, its individual noise intensity is set back to the initial value Sq. 

Figure lb shows the result of computer simulations of trunk trail formation. The related 
dendritic trail system of Pheidole milicida, a harvesting ant of the southwestern U.S. deserts, 
is displayed in Figure la. In our simulation, a nest is assumed in the middle of a triangular 
lattice of size 100 x 100 with periodic boundary conditions. Initially, there are no chemical 
markings on the lattice. At time t = 0, a number Nq of walkers start from the nest with 
a random direction, leaving markings of chemical 0. If a walker disovers a food source by 
chance, it begins to produce chemical 1. Should such a walker find its way back to the nest, it 
activates an additional number of walkers, the recruits, to move out. The maximum number 
of walkers in the simulation is limited to N^nax, which denotes the population size. 



Sa{t) = {1 - k^)[eo + re{t - t^)]^ 




(14) 
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For the food sources, an extended food distribution at the top and bottom hne of the 
lattice is assumed [^. These sources could be exhausted by the visiting walkers, but the 
accidental discovery of new ones in the neighborhood results in a branching of the main trails 
in the vicinity of the food sources and eventually leads to the dendritic structures. The trail 
system observed in Figure lb remains unchanged in its major parts as has been reported 
also in the biological observations of trunk trail formation by ants [^. Nevertheless, some 



minor trails in the vicinity of the food sources slightly shift in the course of time due to 
fluctuations. 

IV. HUMAN TRAIL FORMATION 



Trail formation by pedestrians has been investigated only very recently [Q. It can be 
interpreted as a complex interplay between pedestrian motion, human orientation, and en- 
vironmental changes: On the one hand, pedestrians tend to take the shortest way to their 
destination. On the other hand, they avoid to walk on bumpy ground, since this is uncom- 
fortable. Therefore, they prefer to use existing trails, but they build a new shortcut, if the 
relative detour would be too large. In the latter case they generate a new trail, since foot- 
prints clear some vegetation. Examples of the resulting trail systems can be found in green 
areas, like public parks (cf. Fig. 3). 

Empirical studies have shown that pedestrian motion can be surprisingly well described 



by the social force model sketched in Section II [qT,29|. In particular, it has been demon- 



strated that this model allows a realistic simulation of various observed self-organization 



phenomena in pedestrian crowds p8|-|33|. This includes the emergence of collective patterns 



of motion, e.g. lanes of uniform walking direction p0| , |33| or roundabout traffic at intersec 



tions pTH33|] 



In this section, however, we want to model the evolution of human trail patterns. We will 
assume that the pedestrians behave 'reasonably' and, as before, we will restrict our model 
to the most important factors. It is obvious that pedestrians are able to show a much more 
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complicated behavior than described here. 

Since the equation (||) of motion can be also applied to pedestrians, we have to specify 
how moving pedestrians change their environment by leaving footprints, now. This time, we 
do not have to distinguish different kinds of markings. Thus, we will need only one ground 
potential G{r,t), and the subscript k can be omitted. The value of G is a measure of the 
comfort of walking. (Therefore, it can considerably depend on the weather conditions, which 
is not discussed here any further.) 

For the strength Qair,t) of the markings produced by footprints at place r we assume 

G(r,t) 



Q«(r,t) = /(r) 



1 - 



(15) 



Cfmax (^) 

where J(r) is the location-dependent intensity of clearing vegetation. The saturation term 
[1 — G{r, t)/G.maxi'f)] results from the fact that the clarity of a trail is limited to a maximum 
value Gmax{r). 

On a plain, homogeneous ground without any trails, the desired direction Cq, of a pedes- 
trian a at place r is given by the direction e* of the next destination <i„, i.e. 

e^ir,v,t) = e;«r) = j° ~ = Vf/,(r) (16) 

1 1 ""a 1 1 

with the destination potential 

U^{r) = -\\d^-r\\. (17) 

However, the perception of already existing trails will have an attractive effect f^j.{r,t) 
on the walker, which will again be defined by the gradient of the trail potential Vtr{r,t), 
specified later on: 

f,,{r,t) = VVUr,t). (18) 

Since the potentials U and Vtr influence the pedestrian at the same time, it seems reasonable 
to introduce an orientation relation similar to (|T3|), by taking the sum of both potentials: 
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ea{r,v,t) 




(19) 



Here, N'{r,t) = ||V[f/a(T') + V^r.(r,t)]|| serves as normalization factor. By relation ([l9|) we 
reach that the vector ea{ra,t) points into a direction which is a compromise between the 
shortness of the direct way to the destination and the comfort of using an existing trail. 

Finally, we need to specify the trail potential Vtr for pedestrians. Obviously a trail must 
be recognized by the walkers and near enough in order to be used. Whereas the ground 
potential G{r,t) describes the existence of a trail segment at position r, the trail potential 
Vtr{fa,t) reflects the attractiveness of a trail from the actual position ra{t) of the walker. 
Since this will decrease with the distance ||r — r^H, we have applied the relation 



where a{ra) characterizes the sight, i.e. the range of visibility. In analogy to (p!0|), this 
formula could be easily generalized to include conceivable effects of a pedestrian's angle of 
sight. However, we will not do this here, since we would have to calculate different trail 
potentials for all walkers a, then. This would make the model much more complicated. 

The simulation results of the above described trail formation model are in good agree- 
ment with empirical observations, as can be seen by comparison with photographs. Our 
multi agent simulations begin with plain, homogeneous ground. All pedestrians have their 
own destinations and entry points (like shops, houses, underground stations, or parking 
lots), from which they start at a randomly chosen time. In Figure 4 the entry points and 
destinations are distributed over the small ends of the ground, while in Figure 5 (Fig. 7) 
pedestrians move between all possible pairs of three (four) fixed places. 

At the beginning, pedestrians take the direct ways to their respective destinations. How- 
ever, after some time pedestrians begin to use already existing trails, since this is more 
comfortable than to clear new ways. The frequency of usage decides which trails are rein- 
forced and which ones vanish in the course of time. If the attractiveness of the forming trails 




(20) 
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is large, the final trail system is a minimal way system (which is the shortest way system that 
connects all entry points and destinations). However, because of the pedestrians' dislike of 
taking detours the evolution of the trail system normally stops before this state is reached. 
In other words, a so-called minimal detour system develops if the model parameters are 
chosen realistically (cf. Fig. 5). The resulting trails can considerably differ from the direct 
ways which the pedestrians would use if these were equally comfortable. 



A. Scaling to dimensionless equations 

The use of existing trails depends on the visibility, as given by Eq. (^). Assuming that 
the sight parameter a is approximately space- independent, an additional simplification of 
the equations of trail formation can be reached by introducing dimensionless variables 

x = -, (21) 
a 

G'{x,t) =aG{ax,TT), (23) 

K;(...)^/<^Ve--I.GV'.r). (24, 

U',ix) = -\\dja-x\\, (25) 

etc. Neglecting fluctuations in equation (^ for the moment, this implies the following scaled 
equations: 

— = v^{x^, t) ^ e„(a;„, r) (26) 

CLT a 

for pedestrian motion, 

V[K{x) + V,[.ix,T)] 



e^{x,T) 



\V[U:,{x) + V,[{x,t)]\\ 
for human orientation, and 

G'(x,t) 



(27) 



dT 



1 



G^^^ (x} 



EI(ax)T(ax) ^, , 
^ ^5(x-a;,(r)) (28) 
a 
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for environmental changes. Therefore, we find the surprising result that the dynamics of 
trail formation is (apart from the influence of the number and places of entering and leaving 
pedestrians) already determined by two local parameters k and A instead of four, namely 
the products 

, , I(ax)T(ax) 
k{x) = ^ — '-^ — '- (29) 



and 



= (30) 



Herein, denotes the mean value of the desired velocities v^. 



B. Macroscopic formulation of trail formation 

From the above 'microscopic' model of trail formation we will now derive the related 
'macroscopic' equations. For this purpose we need to distinguish different subpopulations 
a of individuals a. By a(r) we denote the time-dependent set of individuals a who have 
started from the same entry point p„ with the same destination da- Therefore, the different 
sets a correspond to the possible (directed) combinations between existing entry points and 
destinations. 

Next, we define the density Pa{x,T) of individuals of subpopulation a at place x by 

Pa{x,T) = ^ 6{x-Xa{T)). (31) 

aGa(r) 

Note that a spatial smoothing of the density is reached by a discretization of space, which 
is needed for a numerical implementation of the model. For example, if the discrete places 
Xi represent quadratic domains 

A{xi) = {x : \\x - Xi\\^ < L} (32) 

with an area |^| = L^, the corresponding density is 
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Pa{xi,T) = J (fx ^ 5{x-Xo,{t)). (33) 

However, for reasons of simplicity we will treat the continuous case. 
The quantity 

= H j d^^5{x-x^{T)) (34) 

oSa{T) 

describes the number of pedestrians of subpopulation a, who are walking on the ground at 
time r. It changes by pedestrians entering the system at the entry point with a rate 
R^{p^, t) and leaving it at the destinations da with a rate R~{da, t). 
Due to the time-dependence of the sets a(r), we will need the set 

an(r) = a(r + A) n a(r) (35) 
of pedestrians remaining in the system, the set 

a+(r) = a(r + A) \an(r) (36) 

of entering pedestrians, and the set 

a-{T) = a{T)\an{T) (37) 

of leaving pedestrians, for which the following relations hold: 

a+(r) n a_(r) = 0, (38) 
an(r) Ua+(r) = a(r + A), (39) 
an(r) U a_(r) = a(r) . (40) 



Therefore, equation (^) implies 



dpj^ = hm 1 
dr A^o A 



5{x-x^{t + A)) 

aga(r+A) 
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= l!^o ^ ^[^i^ - ^c.{r + A)) 

aean (r) 

-5{x - Xa{r)) 

aea+(r) 



a&a- (r) 



Taking into account 



lim —5(x — xJr + A)) = lim —5(x — Xair)) 

d 

+— 5(a; - Xair)) , 



which follows by Taylor expansion, we obtain 



dpa{x,T) 
dr 



q;6o(t+A) 

-5{x - Xair)) 



+ ^™ol ^ S{x-Xa{T)) 

aea+(T) 



aGa_ (t) 



With 



S{X — Xa{T + A)) — 6{X — Xair)) 



lim — ^ 

A^O A 



and the relations 



aea+ (r) 



^a(iK,T) =hm^ J2 S{x-Xa{r)) 



aEa- (r) 
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for the rates of pedestrians joining and leaving subpopulation a, we finally arrive at 

dpaix,T) 



dr 



»€ci(t) 

+R+ix,T)-R-{x,T), 



(47) 



where R^ix, r) is zero away from the entry point p^, and the same holds for R^ix, r) away 
from the destination da- 

Now, we define the average velocity V a by 



Va{x,T) = v'^{Xc.,t)5{x - Xa{T)) 

Pa{X,T) ^ 



(48) 



This gives us the desired continuity equation 

dpa{x,T 



dr 



+ V ■[pa{x,T)Va{x,T)\ 

= RUx,r)-R-{x,T) 



(49) 



describing pedestrian motion. Fluctuation effects can be taken into account by the additional 

diffusion terms 



Y,^-[Dab{{Pc})Vp,{x,r)] 



(50) 

on the right-hand side of (|49|) [^^. This causes the trails to become somewhat broader. 
Next, we rewrite equation (^) for environmental changes in the form 

dG'ix, t) 



dr 



mx)-G'{x,T)] 



+ 



1 



G'{x,t) 



YK{x)pa{x,t) 



(51) 



Finally, the orientation relation becomes 

V[Ua{x) + VlXx,T)] 



e„{x,T 



\V[Ua{x) + V/,{x,T) 



(52) 



with 



Ua{x) = -\\da/cr - X\ 



(53) 
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Therefore, the average velocity is given by 

Va{x, r) ^ V !^^^^e„(x, t)5{x - x^{t)) 

^^—^^^^^ea{x,T) = \{x)ea{x,T), (54) 
cr 

where is again the average desired pedestrian velocity. In the case of frequent interac- 
tions (avoidance maneuvers) of pedestrians, must be replaced by suitable monotonically 
decreasing functions K({pc}) of the densities pc PH| . pT] , |^ . Moreover, fluctuation effects will 
be stronger, leading to greater diffusion functions Dab{{pc\) and broader trails. 

Summarizing our results, we have found a macroscopic formulation of trail formation 
which is given by equations (^Of) through (^) with (p4|). Apart from possible analytical 



investigations, it allows us to determine implicit equations for the stationary solution, if the 
rates R^{x^ r) and Ra{x, t) are time-independent. Setting the temporal derivatives to zero, 
we find the relations 

^ ^ l + Ea'^NPa(^)/GU.(^) ^ ^ 

and 

V ■ [p,{x)Va{x)] = R+{x) - R-{x) . (56) 
Together with equations (|2^) and (^) through (jS^) , relations ( ^5|) and (|56D allow to calculate 



the finally evolving trail system. Again, we see that the resulting state depends on the two 
parameters A and k. In addition, it is determined by the respective boundary conditions, 
i.e. the configuration and frequency of usage of the entry point-destination pairs, which are 
characterized by the concrete form of the entering rates R'^{x) and leaving rates R~{x). 

The advantage of applying the macroscopic equations is that the finally evolving trail 
system can be calculated much more efficiently, since considerably less time is required for 
computing: The numerical solution can now be obtained by means of a simple iterative 
method which is comparable to the self-consistent field technique. Examples are shown in 
Figure 6 for different values of n. As expected, the results agree with the ones of the related 
microsimulations, which are depicted in Figure 5. 
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We showed that the active walker concept is suitable for modeling and understanding 
trail formation by pedestrians and animals. Our model turned out to be in good agreement 
with observations. It included an equation of motion of the walkers, an equation describing 
environmental changes by the markings which they leave and their decay, a relation reflecting 
the attractiveness of already existing trails, and an equation delineating their influence on 
orientation. Whereas frequently used trails are reinforced, rarely chosen trails are vanishing 
in the course of time. This causes a tendency of trail bundling, which can be interpreted 
as an agglomeration phenomenon. However, the evolving patterns are not localized since 
the active walkers intend to reach certain destinations, starting from their respective entry 
points. 

The structure of the resulting trail system can considerably vary with the species. This 
depends decisively on the main effect which counteracts the trail attraction. Whereas our 
model ants find their destinations (the food sources) by chance, pedestrians can directly 
orient towards their destinations, so that fluctuations are no necessary model component in 
this case. Thus, for certain ant species a dendritic trail system is found, the detailed form 
of which depends on random events, i.e. the concrete history of its evolution. Pedestrians, 
however, produce a minimal detour system, i.e. an optimal compromise between a direct 
way system and a minimal way system. 

As a consequence, we could derive a macroscopic model for the trail formation by pedes- 
trians, but not for ants. It implied a self-consistent field method for a very efficient calculation 
of the finally evolving trail system. This is determined by the location of the entry points 
and destinations (e.g. houses, shops, or parking lots) and the rates of choosing the possible 
connections between them. Apart from this it depends on two parameters only, which was 
demonstrated by scaling to dimensionless equations. These are related to the trail attrac- 
tiveness and the average velocity of motion. 
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In order to demonstrate that the evolution of trail systems can be understood as a typical 
self-organization phenomenon, our model has made a number of simplifying assumptions 
about the agents. 

In the example of trail formation by certain ant species, the major difference to biology 
is, that the active walkers used in the simulations have far less complex capabilities than 
the biological creatures. They almost behave like physical particles which respond to local 
forces in a quite simple manner, without "implicit and explicit intelligence" [^]. Compared 
to the complex 'individual-based' models in ecology [H, the active walker model proposed 
here provides a very simple but efficient tool to simulate a specific structure only with a few 
adjustable parameters. 

With respect to the formation of trunk trails, our model indicates that these patterns 
can be obtained also under the restrictions, that (i) no visual navigation and internal storage 
of information is provided, (ii) in the beginning, no chemical signposts exist which lead the 
ants to the food sources and afterwards back to the nest. Rather, the formation of trail 
systems can be described as a process of self-organization. Based on the interactions of the 
active walkers on a local or 'microscopic' level, the emergence of a global or 'macroscopic' 
structure occurs. The basic interaction between the active walkers can be considered as 
indirect communication mediated by an external storage medium [^,^ . This is a collective 
process in which all active walkers are involved. The information which an active walker 
produces in terms of chemical markings affects the behaviors of the others. It can be amplified 
during the evolution process or disappear again, thus leading to a correlation between the 
information generated and to the self-organization of the walkers on a spatial level. 
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Computer simulations of our pedestrian trail formation model will be a valuable tool for 
designing convenient way systems (cf. Fig. 7). For planning purposes the model parameters 
A and k must be specified in a realistic way. Then, one needs to simulate the exptected flows 
of pedestrians that enter the considered system at certain entry points with the intention to 
reach certain destinations. Already existing ways can be taken into account by the function 
Gq{x). According to our model, a trail system will evolve which minimizes overall detours 
and thereby provides an optimal compromise between a direct and a minimal way system. 
It is expected that the corresponding ways meet the pedestrian requirements best: They 
will most likely be accepted and actually used, since they take into account the route choice 
habits of pedestrians. For the simulation of realistic situations, the results can serve as 
planning guidelines for architects, landscape gardeners, and urban planners. 



C. Current research directions 

Besides of possible applications, our present research focusses on two questions: (i) How 
must our trail formation model be specified in order to be applicable to trail formation by 
hoofed animals or mice |]0^JH3[] ? (ii) Can our model be generalized in a way that allows 
to understand human decision making, in particular processes of finding suitable compro- 
mises? Interestingly enough, one says that someone "follows in somebody's footsteps" or 
that someone "treads new paths". Therefore, a related theory for a more abstract space 
(which represents the set of behavioral alternatives) may describe the evolution of social 



norms and conventions 11 . 
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FIGURES 




FIG. 1. (a) Dendritic trunk trail system of the ant species Pheidole militicida (after [^). (b) 
Simulation result of trunk trail formation by active walkers. The result is in good agreement with 
the empirical findings. 
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FIG. 2. Illustration of the model quantities characterizing the ant-like active walker a. The 
arrow represents the body of the ant and its orientation cua with respect to the axes x and y of 
the co-ordinate system. It ends at the point which corresponds to the front of the ant's head, 
where the antennae start. These have the length Ar and include an angle 2ip of perception. 
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FIG. 3. Between the straight, paved ways on the university campus in Stuttgart- Vaihingen a 
trail system has evolved (center of the picture). Two types of nodes are observed: Intersections of 
two trails running in a straight line and junctions of two trails which smoothly merge into one trail 



|3|,|8|. 
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FIG. 4. When pedestrians leave footprints on the ground, trails will develop, and only parts of 
the ground are used for walking (in contrast to paved areas) . The similarity between the simulation 
result (left) and the trail system on the university campus of Brasilia (right, reproduction by kind 



permission of Klaus Humpert) is obvious |l34|28 |. 
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FIG. 5. The structure of the emerging trail system (hght grey) essentially depends on the 
attractiveness of the trails (i.e. on the parameter k = IT /a). If trail attractiveness is small, a 
direct way system develops (left), if it is large, a minimal way system is formed, otherwise a 
minimal detour system will result (middle) which looks similar to the trail system in the center of 
Fig. 0. The grey scale allows to reconstruct the temporal evolution of the trail system before its 



final state was reached [28 




YY 



FIG. 6. Stationary solution of the macroscopic trail formation model, obtained by an itera- 
tive self-consistent field method. Since the boundary conditions were chosen as in Figure |5|, the 
results in dependence of the parameter n are almost identical with those of the corresponding 
microsimulations 
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FIG. 7. Comparison of different types of way systems between four places: The direct way 
system (which is represented by the black lines) provides the shortest connections between all 
entry points and destinations, but it covers much space. In real situations, pedestrians will produce 
a 'minimal detour system' as the best compromise between a direct way system and a minimal 
way system (which is the shortest way system that connects all entry points and destinations) 
|28|. The illustration shows a simulation result which could serve as a planning guideline. Its 
asymmetry is caused by differences in the frequency of trail usage. (Note that the above figure, 
in contrast to Figs. ^ and ^, does not display the ground potential, but the trail potential. The 
latter appears considerably broader, since it takes into account the range of visibility of the trails. 
Arrows represent the positions and walking directions of pedestrians. Therefore, they indicate the 
actually taken ways.) 



